CRISPR-broad: combined design of multi-targeting gRNAs and broad, multiplex target finding.

In CRISPR-Cas and related nuclease-mediated genome editing, target recognition is based on guide RNAs (gRNAs) that are complementary to selected DNA regions. While single site targeting is fundamental for localized genome editing, targeting to expanded and multiple chromosome elements is desirable for various biological applications such as genome mapping and epigenome editing that make use of different fusion proteins with enzymatically dead Cas9. The current gRNA design tools are not suitable for this task, as these are optimized for defining single gRNAs for unique loci. Here, we introduce CRISPR-broad, a standalone, open-source application that defines gRNAs with multiple but specific targets in large continuous or spread regions of the genome, as defined by the user. This ability to identify multi-targeting gRNAs and corresponding multiple targetable regions in genomes is based on a novel aggregate gRNA scoring derived from on-target windows and off-target sites. Applying the new tool to the genomes of two model species, C. elegans and H. sapiens, we verified its efficiency in determining multi-targeting gRNAs and ranking potential target regions optimized for broad targeting. Further, we demonstrated the general usability of CRISPR-broad by cellular mapping of a large human genome element using dCas9 fused to green fluorescent protein.

Microbe-induced drought tolerance by ABA-mediated root architecture and epigenetic reprogramming.

The use of beneficial microbes to mitigate drought stress tolerance of plants is of great potential albeit little understood. We show here that a root endophytic desert bacterium, Pseudomonas argentinensis strain SA190, enhances drought stress tolerance in Arabidopsis. Transcriptome and genetic analysis demonstrate that SA190-induced root morphogenesis and gene expression is mediated via the plant abscisic acid (ABA) pathway. Moreover, we demonstrate that SA190 primes the promoters of target genes in an epigenetic ABA-dependent manner. Application of SA190 priming on crops is demonstrated for alfalfa, showing enhanced performance under drought conditions. In summary, a single beneficial root bacterial strain can help plants to resist drought conditions.

PhaeoEpiView: an epigenome browser of the newly assembled genome of the model diatom Phaeodactylum tricornutum.

Recent advances in DNA sequencing technologies particularly long-read sequencing, greatly improved genomes assembly. However, this has created discrepancies between published annotations and epigenome tracks, which have not been updated to keep pace with the new assemblies. Here, we used the latest improved telomere-to-telomere assembly of the model pennate diatom Phaeodactylum tricornutum to lift over the gene models from Phatr3, a previously annotated reference genome. We used the lifted genes annotation and newly published transposable elements to map the epigenome landscape, namely DNA methylation and post-translational modifications of histones. This provides the community with PhaeoEpiView, a browser that allows the visualization of epigenome data and transcripts on an updated and contiguous reference genome, to better understand the biological significance of the mapped data. We updated previously published histone marks with a more accurate peak calling using mono instead of poly(clonal) antibodies and deeper sequencing. PhaeoEpiView ( https://PhaeoEpiView.univ-nantes.fr ) will be continuously updated with the newly published epigenomic data, making it the largest and richest epigenome browser of any stramenopile. In the upcoming era of molecular environmental studies, where epigenetics plays a significant role, we anticipate that PhaeoEpiView will become a widely used tool.

The Promising Epigenetic Regulators for Refractory Epilepsy: An Adventurous Road Ahead.

The attribution of seizure freedom is yet to be achieved for patients suffering from refractory epilepsy, e.g. Dravet Syndrome (DS). The confined ability of mono-chemical entity-based antiseizure drugs (ASDs) to act directly at genomic level is one of the factors, combined with undetermined seizure triggers lead to recurrent seizure (RS) in DS, abominably affecting the sub-genomic architecture of neural cells. Thus, the RS and ASD appear to be responsible for the spectrum of exorbitant clinical pathology. The RS distresses the 5-HT-serotonin pathway, hypomethylates genes of CNS, and modulates the microRNA (miRNA)/long non-coding RNA (lncRNA), eventually leading to frozen molecular alterations. These changes shall be reverted by compatible epigenetic regulators (EGR) like, miRNA and lncRNA from Breast milk (BML) and Bacopa monnieri (BMI). The absence of studious seizure in SCN1A mutation-positive babies for the first 6 months raises the possibility that the consequences of mutation in SCN1A are subsidized by EGRs from BML. EGR-dependent-modifier gene effect is likely imposed by the other members of the SCN family. Therefore, we advocate that miRNA/lncRNA from BML and bacosides/miRNA from BMI buffer the effect of SCN1A mutation by sustainably maintaining modifier gene effect in the aberrant neurons. The presence of miRNA-155-5p, -30b-5p, and -30c-5p family in BML and miR857, miR168, miR156, and miR158 in BMI target at regulating SCN family and CLCN5 as visualized by Cystoscope. Thus, we envisage that the possible effects of EGR might include (a) upregulating the haploinsufficient SCN1A strand, (b) down-regulating seizure-elevated miRNA, (c) suppressing the seizure-induced methyltransferases, and (d) enhancing the GluN2A subunit of NMDA receptor to improve cognition. The potential of these EGRs from BML and BML is to further experimentally strengthen, long-haul step forward in molecular therapeutics.

Exploring the biological behavior of Heat shock protein (HSPs) for understanding the Anti-ischemic stroke in humans.

BACKGROUND: An ischemic stroke can be caused by thrombosis and ischemia, which is a major public health problem around the world, resulting in severe disability and a high death rate. The goal of this work is to examine and target various heat shock proteins (HSPs) via their interacting partners, which may have an anti-ischemic stroke impact. METHODS: Various heat shock proteins are identified and used for construction of PPI network through STRING webserver. Networks are analysed and visualized using the cytoscape for checking the protein-protein interactions. Along with this, multiple cytoscape based modules are integrated for the analysis of results, and Gene Ontology results are analysed using GOView. RESULTS: The core PPI network was revealed with 129 nodes and 1174 edges. Through Gene ontology (GO) and KEGG enrichment analysis the promising function of HSPs in two important signaling pathways were mainly recorded, representing the HSPs are necessary for repair and activations of brain cells during ischemic stroke. In addition, the study is revelation for targeting multiple HSPs via their interacting partners, which can provide anti-ischemic stroke effect. CONCLUSION: Overall, this finding provides a network-based framework for future research on HSP as therapeutic molecules for anti-ischemic stroke related applications.

Evaluating the structural and immune mechanism of Interleukin-6 for the investigation of goat milk peptides as potential treatments for COVID-19.

The function of Immune control, haematopoiesis, and inflammation all depend on the cytokine Interleukin 6 (IL-6), and higher expression of IL-6 is seen in COVID-19 and other diseases. The immune protein IL-6 activation is dependent on binding interactions with IL-6Rα, mIL-6R, and sIL-6R for its cellular function. Termination of these reaction could benefit for controlling the over-expression in COVID-19 patients and that may arise as inhibitors for controlling COVID-19. Traditionally, the goat milk has been prescribed as medicine in ayurvedic practice and through this work, we have explored the benefits of peptides from goat milk as IL-6 inhibitors, and it have the potential of inhibiting the over expression of IL-6 and control the COVID-19 disease. Computational experiments have shown that goat peptides had strong interactions with IL-6, with higher scoring profiles and energy efficiency ranging from -6.00 kcal/mol to -9.00 kcal/mol in docking score and -39.00 kcal/mol in binding energy. Especially the YLGYLEQLLR, VLVLDTDYK and AMKPWIQPK peptides from goat milk holds better scoring and shows strong interactions were identified as the most potential IL-6 inhibitor candidates in this study. Peptides from Goat proteins, which are capable of binding to the IL-6 receptor with strong binding conformations, have no negative effects on other immune system proteins.

Information theoretic perspective on genome clustering.

Shannon's information theoretic perspective of communication helps one to understand the storage and processing of information in one-dimensional sequences. An information theoretic analysis of 937 available completely sequenced prokaryotic genomes and 238 eukaryotic chromosomes is presented. Information content (Id) values were used to cluster these chromosomes. Chargaff's second parity rule i.e compositional self-complementarity, an empirical fact is observed in all the genomes, except for the proteobacteria Candidatus Hodgkinia cicadicola. High information content, arising out of biased base composition in all the 14 chromosomes of Plasmodium falciparum is found among two other genomes of prokaryotes viz. Buchnera aphidicola str. Cc (Cinara cedri) and Candidatus Carsonella ruddii PV. Despite size and compositional variations, both prokaryotic and eukaryotic genomes do not deviate significantly from an equiprobable and random situation. Eukaryotic chromosomes of an organism tend to have similar informational restraints as seen when a simple distance based method is used to cluster them. In eukaryotes, in certain cases, Id values are also similar for the two arms (p and q arm) of the chromosomes. The results of this current study confirm that the information content can provide insights into the clustering of genomes and the evolution of messaging strategies of the genomes. An efficient and robust Perl CGI standalone tool is created based on this information theory algorithm for the analysis of the whole genomes and is made available at https://github.com/AlagurajVeluchamy/InformationTheory.

Genome wide natural variation of H3K27me3 selectively marks genes predicted to be important for cell differentiation in Phaeodactylum tricornutum.

In multicellular organisms, Polycomb Repressive Complex2 (PRC2) is known to deposit tri-methylation of lysine 27 of histone H3 (H3K27me3) to establish and maintain gene silencing, critical for developmentally regulated processes. The PRC2 complex is absent in both widely studied model yeasts, which initially suggested that PRC2 arose with the emergence of multicellularity. However, its discovery in several unicellular species including microalgae questions its role in unicellular eukaryotes. Here, we use Phaeodactylum tricornutum enhancer of zeste E(z) knockouts and show that P. tricornutum E(z) is responsible for di- and tri-methylation of lysine 27 of histone H3. H3K27me3 depletion abolishes cell morphology in P. tricornutum providing evidence for its role in cell differentiation. Genome-wide profiling of H3K27me3 in fusiform and triradiate cells further revealed genes that may specify cell identity. These results suggest a role for PRC2 and its associated mark in cell differentiation in unicellular species, and highlight their ancestral function in a broader evolutionary context than currently is appreciated.

Alternative splicing and allosteric regulation modulate the chromatin binding of UHRF1.

UHRF1 is an important epigenetic regulator associated with apoptosis and tumour development. It is a multidomain protein that integrates readout of different histone modification states and DNA methylation with enzymatic histone ubiquitylation activity. Emerging evidence indicates that the chromatin-binding and enzymatic modules of UHRF1 do not act in isolation but interplay in a coordinated and regulated manner. Here, we compared two splicing variants (V1, V2) of murine UHRF1 (mUHRF1) with human UHRF1 (hUHRF1). We show that insertion of nine amino acids in a linker region connecting the different TTD and PHD histone modification-binding domains causes distinct H3K9me3-binding behaviour of mUHRF1 V1. Structural analysis suggests that in mUHRF1 V1, in contrast to V2 and hUHRF1, the linker is anchored in a surface groove of the TTD domain, resulting in creation of a coupled TTD-PHD module. This establishes multivalent, synergistic H3-tail binding causing distinct cellular localization and enhanced H3K9me3-nucleosome ubiquitylation activity. In contrast to hUHRF1, H3K9me3-binding of the murine proteins is not allosterically regulated by phosphatidylinositol 5-phosphate that interacts with a separate less-conserved polybasic linker region of the protein. Our results highlight the importance of flexible linkers in regulating multidomain chromatin binding proteins and point to divergent evolution of their regulation.

Genome-enabled phylogenetic and functional reconstruction of an araphid pennate diatom Plagiostriata sp. CCMP470, previously assigned as a radial centric diatom, and its bacterial commensal.

Diatoms are an ecologically fundamental and highly diverse group of algae, dominating marine primary production in both open-water and coastal communities. The diatoms include both centric species, which may have radial or polar symmetry, and the pennates, which include raphid and araphid species and arose within the centric lineage. Here, we use combined microscopic and molecular information to reclassify a diatom strain CCMP470, previously annotated as a radial centric species related to Leptocylindrus danicus, as an araphid pennate species in the staurosiroid lineage, within the genus Plagiostriata. CCMP470 shares key ultrastructural features with Plagiostriata taxa, such as the presence of a sternum with parallel striae, and the presence of a highly reduced labiate process on its valve; and this evolutionary position is robustly supported by multigene phylogenetic analysis. We additionally present a draft genome of CCMP470, which is the first genome available for a staurosiroid lineage. 270 Pfams (19%) found in the CCMP470 genome are not known in other diatom genomes, which otherwise does not hold big novelties compared to genomes of non-staurosiroid diatoms. Notably, our DNA library contains the genome of a bacterium within the Rhodobacterales, an alpha-proteobacterial lineage known frequently to associate with algae. We demonstrate the presence of commensal alpha-proteobacterial sequences in other published algal genome and transcriptome datasets, which may indicate widespread and persistent co-occurrence.

GCN5 modulates salicylic acid homeostasis by regulating H3K14ac levels at the 5' and 3' ends of its target genes.

The modification of histones by acetyl groups has a key role in the regulation of chromatin structure and transcription. The Arabidopsis thaliana histone acetyltransferase GCN5 regulates histone modifications as part of the Spt-Ada-Gcn5 Acetyltransferase (SAGA) transcriptional coactivator complex. GCN5 was previously shown to acetylate lysine 14 of histone 3 (H3K14ac) in the promoter regions of its target genes even though GCN5 binding did not systematically correlate with gene activation. Here, we explored the mechanism through which GCN5 controls transcription. First, we fine-mapped its GCN5 binding sites genome-wide and then used several global methodologies (ATAC-seq, ChIP-seq and RNA-seq) to assess the effect of GCN5 loss-of-function on the expression and epigenetic regulation of its target genes. These analyses provided evidence that GCN5 has a dual role in the regulation of H3K14ac levels in their 5' and 3' ends of its target genes. While the gcn5 mutation led to a genome-wide decrease of H3K14ac in the 5' end of the GCN5 down-regulated targets, it also led to an increase of H3K14ac in the 3' ends of GCN5 up-regulated targets. Furthermore, genome-wide changes in H3K14ac levels in the gcn5 mutant correlated with changes in H3K9ac at both 5' and 3' ends, providing evidence for a molecular link between the depositions of these two histone modifications. To understand the biological relevance of these regulations, we showed that GCN5 participates in the responses to biotic stress by repressing salicylic acid (SA) accumulation and SA-mediated immunity, highlighting the role of this protein in the regulation of the crosstalk between diverse developmental and stress-responsive physiological programs. Hence, our results demonstrate that GCN5, through the modulation of H3K14ac levels on its targets, controls the balance between biotic and abiotic stress responses and is a master regulator of plant-environmental interactions.

Wheat chromatin architecture is organized in genome territories and transcription factories.

BACKGROUND: Polyploidy is ubiquitous in eukaryotic plant and fungal lineages, and it leads to the co-existence of several copies of similar or related genomes in one nucleus. In plants, polyploidy is considered a major factor in successful domestication. However, polyploidy challenges chromosome folding architecture in the nucleus to establish functional structures. RESULTS: We examine the hexaploid wheat nuclear architecture by integrating RNA-seq, ChIP-seq, ATAC-seq, Hi-C, and Hi-ChIP data. Our results highlight the presence of three levels of large-scale spatial organization: the arrangement into genome territories, the diametrical separation between facultative and constitutive heterochromatin, and the organization of RNA polymerase II around transcription factories. We demonstrate the micro-compartmentalization of transcriptionally active genes determined by physical interactions between genes with specific euchromatic histone modifications. Both intra- and interchromosomal RNA polymerase-associated contacts involve multiple genes displaying similar expression levels. CONCLUSIONS: Our results provide new insights into the physical chromosome organization of a polyploid genome, as well as on the relationship between epigenetic marks and chromosome conformation to determine a 3D spatial organization of gene expression, a key factor governing gene transcription in polyploids.

R-Loop Mediated trans Action of the APOLO Long Noncoding RNA.

In eukaryotes, three-dimensional genome organization is critical for transcriptional regulation of gene expression. Long noncoding RNAs (lncRNAs) can modulate chromatin conformation of spatially related genomic locations within the nucleus. Here, we show that the lncRNA APOLO (AUXIN-REGULATED PROMOTER LOOP) recognizes multiple distant independent loci in the Arabidopsis thaliana genome. We found that APOLO targets are not spatially associated in the nucleus and that APOLO recognizes its targets by short sequence complementarity and the formation of DNA-RNA duplexes (R-loops). The invasion of APOLO to the target DNA decoys the plant Polycomb Repressive Complex 1 component LHP1, modulating local chromatin 3D conformation. APOLO lncRNA coordinates the expression of distal unrelated auxin-responsive genes during lateral root development in Arabidopsis. Hence, R-loop formation and chromatin protein decoy mediate trans action of lncRNAs on distant loci. VIDEO ABSTRACT.

A genomics approach reveals the global genetic polymorphism, structure, and functional diversity of ten accessions of the marine model diatom Phaeodactylum tricornutum.

Diatoms emerged in the Mesozoic period and presently constitute one of the main primary producers in the world's ocean and are of a major economic importance. In the current study, using whole genome sequencing of ten accessions of the model diatom Phaeodactylum tricornutum, sampled at broad geospatial and temporal scales, we draw a comprehensive landscape of the genomic diversity within the species. We describe strong genetic subdivisions of the accessions into four genetic clades (A-D) with constituent populations of each clade possessing a conserved genetic and functional makeup, likely a consequence of the limited dispersal of P. tricornutum in the open ocean. We further suggest dominance of asexual reproduction across all the populations, as implied by high linkage disequilibrium. Finally, we show limited yet compelling signatures of genetic and functional convergence inducing changes in the selection pressure on many genes and metabolic pathways. We propose these findings to have significant implications for understanding the genetic structure of diatom populations in nature and provide a framework to assess the genomic underpinnings of their ecological success and impact on aquatic ecosystems where they play a major role. Our work provides valuable resources for functional genomics and for exploiting the biotechnological potential of this model diatom species.

INDETERMINATE-DOMAIN 4 (IDD4) coordinates immune responses with plant-growth in Arabidopsis thaliana.

INDETERMINATE DOMAIN (IDD)/ BIRD proteins are a highly conserved plant-specific family of transcription factors which play multiple roles in plant development and physiology. Here, we show that mutation in IDD4/IMPERIAL EAGLE increases resistance to the hemi-biotrophic pathogen Pseudomonas syringae, indicating that IDD4 may act as a repressor of basal immune response and PAMP-triggered immunity. Furthermore, the idd4 mutant exhibits enhanced plant-growth indicating IDD4 as suppressor of growth and development. Transcriptome comparison of idd4 mutants and IDD4ox lines aligned to genome-wide IDD4 DNA-binding studies revealed major target genes related to defense and developmental-biological processes. IDD4 is a phospho-protein that interacts and becomes phosphorylated on two conserved sites by the MAP kinase MPK6. DNA-binding studies of IDD4 after flg22 treatment and with IDD4 phosphosite mutants show enhanced binding affinity to ID1 motif-containing promoters and its function as a transcriptional regulator. In contrast to the IDD4-phospho-dead mutant, the IDD4 phospho-mimicking mutant shows altered susceptibility to PstDC3000, salicylic acid levels and transcriptome reprogramming. In summary, we found that IDD4 regulates various hormonal pathways thereby coordinating growth and development with basal immunity.

Thermopriming triggers splicing memory in Arabidopsis.

Abiotic and biotic stresses limit crop productivity. Exposure to a non-lethal stress, referred to as priming, can allow plants to survive subsequent and otherwise lethal conditions; the priming effect persists even after a prolonged stress-free period. However, the molecular mechanisms underlying priming are not fully understood. Here, we investigated the molecular basis of heat-shock memory and the role of priming in Arabidopsis thaliana. Comprehensive analysis of transcriptome-wide changes in gene expression and alternative splicing in primed and non-primed plants revealed that alternative splicing functions as a novel component of heat-shock memory. We show that priming of plants with a non-lethal heat stress results in de-repression of splicing after a second exposure to heat stress. By contrast, non-primed plants showed significant repression of splicing. These observations link 'splicing memory' to the ability of plants to survive subsequent and otherwise lethal heat stress. This newly discovered priming-induced splicing memory may represent a general feature of heat-stress responses in plants and other organisms as many of the key components are conserved among eukaryotes. Furthermore, this finding could facilitate the development of novel approaches to improve plant survival under extreme heat stress.

MAPK-triggered chromatin reprogramming by histone deacetylase in plant innate immunity.

BACKGROUND: Microbial-associated molecular patterns activate several MAP kinases, which are major regulators of the innate immune response in Arabidopsis thaliana that induce large-scale changes in gene expression. Here, we determine whether microbial-associated molecular pattern-triggered gene expression involves modifications at the chromatin level. RESULTS: Histone acetylation and deacetylation are major regulators of microbial-associated molecular pattern-triggered gene expression and implicate the histone deacetylase HD2B in the reprogramming of defence gene expression and innate immunity. The MAP kinase MPK3 directly interacts with and phosphorylates HD2B, thereby regulating the intra-nuclear compartmentalization and function of the histone deacetylase. CONCLUSIONS: By studying a number of gene loci that undergo microbial-associated molecular pattern-dependent activation or repression, our data reveal a mechanistic model for how protein kinase signaling directly impacts chromatin reprogramming in plant defense.

Arabidopsis ATRX Modulates H3.3 Occupancy and Fine-Tunes Gene Expression.

Histones are essential components of the nucleosome, the major chromatin subunit that structures linear DNA molecules and regulates access of other proteins to DNA. Specific histone chaperone complexes control the correct deposition of canonical histones and their variants to modulate nucleosome structure and stability. In this study, we characterize the Arabidopsis thaliana Alpha Thalassemia-mental Retardation X-linked (ATRX) ortholog and show that ATRX is involved in histone H3 deposition. Arabidopsis ATRX mutant alleles are viable, but show developmental defects and reduced fertility. Their combination with mutants of the histone H3.3 chaperone HIRA (Histone Regulator A) results in impaired plant survival, suggesting that HIRA and ATRX function in complementary histone deposition pathways. Indeed, ATRX loss of function alters cellular histone H3.3 pools and in consequence modulates the H3.1/H3.3 balance in the cell. H3.3 levels are affected especially at genes characterized by elevated H3.3 occupancy, including the 45S ribosomal DNA (45S rDNA) loci, where loss of ATRX results in altered expression of specific 45S rDNA sequence variants. At the genome-wide scale, our data indicate that ATRX modifies gene expression concomitantly to H3.3 deposition at a set of genes characterized both by elevated H3.3 occupancy and high expression. Together, our results show that ATRX is involved in H3.3 deposition and emphasize the role of histone chaperones in adjusting genome expression.

The Arabidopsis SWI/SNF protein BAF60 mediates seedling growth control by modulating DNA accessibility.

BACKGROUND: Plant adaptive responses to changing environments involve complex molecular interplays between intrinsic and external signals. Whilst much is known on the signaling components mediating diurnal, light, and temperature controls on plant development, their influence on chromatin-based transcriptional controls remains poorly explored. RESULTS: In this study we show that a SWI/SNF chromatin remodeler subunit, BAF60, represses seedling growth by modulating DNA accessibility of hypocotyl cell size regulatory genes. BAF60 binds nucleosome-free regions of multiple G box-containing genes, opposing in cis the promoting effect of the photomorphogenic and thermomorphogenic regulator Phytochrome Interacting Factor 4 (PIF4) on hypocotyl elongation. Furthermore, BAF60 expression level is regulated in response to light and daily rhythms. CONCLUSIONS: These results unveil a short path between a chromatin remodeler and a signaling component to fine-tune plant morphogenesis in response to environmental conditions.

The quest for epigenetic regulation underlying unisexual flower development in Cucumis melo.

BACKGROUND: Melon (Cucumis melo) is an important vegetable crop from the Cucurbitaceae family and a reference model specie for sex determination, fruit ripening and vascular fluxes studies. Nevertheless, the nature and role of its epigenome in gene expression regulation and more specifically in sex determination remains largely unknown. RESULTS: We have investigated genome wide H3K27me3 and H3K9ac histone modifications and gene expression dynamics, in five melon organs. H3K9ac and H3K27me3 were mainly distributed along gene-rich regions and constrained to gene bodies. H3K9ac was preferentially located at the TSS, whereas H3K27me3 distributed uniformly from TSS to TES. As observed in other species, H3K9ac and H3K27me3 correlated with high and low gene expression levels, respectively. Comparative analyses of unisexual flowers pointed out sex-specific epigenetic states of TFs involved in ethylene response and flower development. Chip-qPCR analysis of laser dissected carpel and stamina primordia, revealed sex-specific histone modification of MADS-box genes. Using sex transition mutants, we demonstrated that the female promoting gene, CmACS11, represses the expression of the male promoting gene CmWIP1 via deposition of H3K27me3. CONCLUSIONS: Our findings reveal the organ-specific landscapes of H3K9ac and H3K27me3 in melon. Our results also provide evidence that the sex determination genes recruit histone modifiers to orchestrate unisexual flower development in monoecious species.